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ABSTRACT 

Context. The Gaia satellite will measure highly accurate absolute parallaxes of hundreds of millions of stars by comparing the 
parallactic displacements in the two fields of view of the optical instrument. The requirements on the stability of the 'basic angle' 
between the two fields are correspondingly strict, and possible variations (on the microarcsec level) are therefore monitored by an 
on-board metrology system. Nevertheless, since even very small periodic variations of the basic angle might cause a global offset of 
the measured parallaxes, it is important to find independent verification methods. 

Aims. We investigate the potential use of Galactic Cepheids as standard candles for verifying the Gaia parallax zero point. 
Methods. We simulate the complete population of Galactic Cepheids and their observations by Gaia. Using the simulated data, 
simultaneous fits are made of the parameters of the period-luminosity relation and a global parallax zero point. 
Results. The total number of Galactic Cepheids is estimated at about 20 000, of which nearly half could be observed by Gaia. In the 
most favourable circumstances, including negligible intrinsic scatter and extinction errors, the determined parallax zero point has an 
uncertainty of 0.2 microarcsec. With more realistic assumptions the uncertainty is several times larger, and the result is \ery sensitive 
to errors in the applied extinction corrections. 

Conclusions. The use of Galactic Cepheids alone will not be sufficient to determine a possible parallax zero-point error to the full 
potential systematic accuracy of Gaia. The global verification of Gaia parallaxes will most likely depend on a combination of many 
different methods, including this one. 

Key words. Stars: variables: Cepheids - Space vehicles: instruments - Parallaxes 



1. Introduction 

The Gaia satellite, due for launch in 2013, will measure the 
trigonometric parallaxes of roughly a billion objects in the 
Gal axy and beyond w ith accuracies reaching 10 juas (microarc- 
sec; iLindegrenI |20 1 Oh . This huge improvement over its prede- 
cessor Hipparcos, in terms of accuracy, limiting magnitude and 
number of objects, will revolutionize many areas of stellar and 
Galactic astrophysics (Ferryman et al. 2001). Moreover, it will 
allow entirely new kinds of investigations that depend on the sta- 
tistical combination of very large data sets. One such example is 
the determination of the distance to the Large Magellanic Cloud 
(LMC). The LMC distance is fundamental for the extragalactic 
distance scale, and current estimates put it at 5 kpc (or 20 Aias 
parallax) with a rel ative uncertainty of 5% (iFreedman et al.l 
l200ltlSchaefe3l2008l) . Gaia should be able to observe some ^ 
10^ stars brighter than 20th magnitude in the LMC, with a stan- 
dard error in the individual parallaxes of about 200 //as or better 
Potentially, therefore, the mean LMC parallax as estimated from 
Gaia data could have an accuracy of 200 A^"'^ - 0.06 juas, equiv- 
alent to a relative error in distance of 0.3%, or 0.006 in distance 
modulus. Needless to say, such a result will be extremely in- 
teresting. However, to achieve the A""^ improvement for large 
N requires (1) that the individual par allax errors are e ffectively 
uncorrelated, which may be the case (iHoll et al.ll20"T0l) . and (2) 
that there is no significant global zero-point error (bias) in the 
measured parallaxes. This and many similar examples show that 
even a bias < 0. 1 //as must be considered significant in compar- 
ison with the potential capabilities of the mission. 



Achieving the desired parallax accuracy requires an exceed- 
ingly stable optical instrument in the Gaia satellite. Even ex- 
tremely small periodic variations in the so-called basic angle be- 
tween the two fields of view could l ead to an undes irable global 
offset of the measured parallaxes (lMignardll201 ll) . One well- 
known cause of such variations is the variable heating from solar 
radiation as the satellite rotates. Although the satellite has been 
carefully designed to minimize these effects, the resulting vari- 
ations in the basic angle cannot be completely eliminated. With 
the use of an on-board laser interferometer they will however be 
continuously measured and taken into account in the instrument 
calibration model. However, it is obviously of great importance 
to be able to verify the resulting parallax zero point by indepen- 
dent, astrophysical means. 

A large number of methods are in principle available for as- 
trophysical verification of the Gaia parallaxes. We may distin- 
guish three main classes of methods: (1) A priori knowledge of 
the parallax, using for example quasars that are so distant that 
their parallaxes for the present purpose can be considered to be 
zero. (2) Distances determined by geometric principles not rely- 
ing on trigonometric parallax. Several of these methods depend 
on a combination of doppler velocity, time and angular measure- 
ments, for example orbital parall axes for spectrosco pic binaries 
with an astrometric orbit (e.g., ITorres et al.l Il997l) . expansion 
parallaxes for supernova rem nants and planetary nebulae (e.g., 
lTrimblelll973"tlLi et al.l l2002'). a geometric variant of the Baade- 
Wesselink method for pulsating stars (e.g., Lang^et al. 20021), 
kinematic distances to globular clusters (e.g., van de Ven et al.l 
120061 : [McLaughlin et al.H2006i) . and rotational parallaxes for ex- 
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ternal galaxies (IOlling|l2007h . (3) Distance ratio methods: the 
ratio of the distances to two or more objects equals the inverse 
ratio of their parallaxes. This equality is violated in the presence 
of a parallax zero point error, which can therefore be derived, 
e.g., from photometric distance ratios established by means of 
standard candles, provided that extinction effects can be mas- 
tered. 

Quasars (belonging to class 1 above) are among the more 
promising candidates. It is estimated that Gaia will observe 
around 500000 quasars with individual parallax uncertainties 
around 250-300 fias d ue to their faintness in the Gaia G band 
dLindegren et al.l l2008'). This would lead to an uncertainty in the 
mean zero-point around 0.4 yuas, with a possible additional bias 
from foreground stars contaminating the sample. 

Distance determinations using various geometric principles 
(class 2 above) will be very important for checking the consis- 
tency and reliability of Gaia parallaxes. However, they are not 
entirely free of model assumptions, and are therefore in many 
respects problematic as a means of verifying the Gaia data. 
Moreover, it is doubtful if the number of objects available and 
the achievable accuracies are high enough for the present pur- 
pose. 

In this paper we concentrate on the distance ratio method, 
based on Classical Cepheids as one of the most reliable stan- 
dard candles. In particular, we focus on the Galactic Cepheids, 
which are fewer than the quasars and extragalactic Cepheids but 
with individually more accurate parallaxes. In this method, to 
avoid a circular argument, we make a simultaneous calibration of 
the period-luminosity (P-L) relation and the parallax zero-point. 
Since Gaia will observe at least ten times the ~ 800 currently 
known Galactic Cepheids, we need to create a synthetic popula- 
tion of Galactic Cepheids with the appropriate properties before 
we can simulate the Gaia observations and make a statistical in- 
vestigation of the expected errors. Since it is difficult to quantify 
how much of the final errors will depend on various modelling 
eiTors (say, from a possible non-linearity of the P-L relation), our 
strategy is to consider first a best-case (i.e., optimistic) scenario, 
where all such effects are negligible, and then investigate how 
sensitive the results are to the various assumptions. 

The method of absolute parallax measurements with Gaia is 
closely related to the physical origin of a possible parallax bias, 
and these aspects of the mission are therefore briefly explained in 
Sect. 121 Our modelling of the Galactic Cepheids and their obser- 
vations by Gaia are described in Sect. |3] The results of the simu- 
lation and parameter fitting experiments are discussed in Sect. HI 
followed by the conclusions in Sect.|5] 

2. Cause of a possible parallax bias in Gaia 

In contrast to ground-based parallaxes, which are always mea- 
sured relative to background objects, astrometric satellites such 
as Hipparcos and Gaia in principle allow to determine absolute 
parallaxes thanks to the large difference in p arallax factor be- 
tween the two widely separate d fields of view (lLindegrenl2005l: 
iLindegren & de Bruiinel2005l) . However, this capability depends 
critically on the short-term (~ few hours) stabil ity of the so- 
called basic ang le between the two fields of view (Ivan LeeuwenI 
120051: lMignardll201lh . The optical instrument of Gaia is de- 
signed to be stable on these time scales to within a few //as, 
and an on-board inteiferometric metrology system, the Basic 
Angle Monitor (BAM), will moreover measure short-term vari- 
ation s of the basic angle with a precision < 1 //as every few min- 
utes (ILindegren et al. 2008). The technical design thus guaran- 
tees that measurement biases related to basic-angle variations 
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Fig. 1. Illustrating the effect of parallax on stellar images in the 
preceding (p) and following (/) field of view of Gaia, when the 
Sun is in two different positions (si and S2) relative to the stars. 
z is the spin axis of Gaia, and the basic angle is marked between 
the two fields of view. The displacement due to parallax is al- 
ways directed towards the Sun, as indicated by the thick aiTows 
at p and /, and is measured by Gaia as projected along the scan- 
ning great circle. See Sect. |2]for further explanation. 



are practically negligible in comparison to the random measure- 
ment errors (> 30 //as per field-of-view crossing). Nevertheless, 
when averaging over many stars we should be concerned about 
systematic effects that are much smaller than the random errors. 

The geometry of the observations with respect to the Sun is 
important both for the determination of parallax and for possi- 
ble thennal^erturbations of the instrument. Indeed, as discussed 
bv lvan Leeuw en (2005), a certain systematic variation of the ba- 
sic angle, depending on the satellite spin phase relative to the 
Sun, has almost the same effect on the measurements as a global 
zero-point shift of the parallaxes. Thus, if such a variation ex- 
ists in the real Gaia instrument, and is not recognized by the 
on-board metrology, then the result will be a global bias of the 
derived parallaxes. The cause of this shift can be understood by 
means of Fig. [T] Stellar parallax causes an apparent shift of the 
star along a great circle towards the Sun. If all stars have a con- 
stant positive parallax, the apparent shifts will be as indicated by 
the short arrows in the figure, depending on the position of the 
Sun relative to the stars. Thanks to its one-dimensional measure- 
ment principle, Gaia is only sensitive to the relative shift along 
the scanning great circle through p and /. Thus, with the Sun 
at Si in the diagram, i.e., closest to the fields of view, the stellar 
parallax shifts on the detectors will be indistinguishable from a 
slight enlargement of the basic angle. With the Sun at S2, furthest 
away from the fields, the stellar shifts will be indistinguishable 
from a slight reduction of the basic angle. Thus, a temporal vari- 
ation of the basic angle, caused for example by the solar heating, 
would mimic a global shift of the parallaxes. 

According to the technical specifications of Gaia, the ex- 
pected parallax bias from this effect is at most a few //as, and 
most of it will be removed in the data processing by means of the 
BAM metrology data. Indeed, the BAM specifications are such 
that, theoretically, the remaining parallax bias should be much 
less than 0.1 //as. Nevertheless, due to the subtle nature of the 
effect, it is important to verify the parallax zero point by some 
independent method. Since the relevant instrumental effects are 
almost completely correlated with a parallax shift, an indepen- 
dent verification must be based on astrophysical considerations. 
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Fig. 2. The upper panel shows the period distribution for the 
Cepheids in the full Berdnikov catalogue (455 Cepheids) as well 
as the volume-complete sample within the assumed complete- 
ness limit of rcosb < 1 kpc (71 Cepheids). Note the larger 
fraction of long-period Cepheids in the full Berdnikov catalogue 
due to selection effects. The line histogram in lower panel shows 
our modelled period distribution consisting of two overlapping 
Gaussians, which ideally should represent the volume-complete 
sample in the upper panel. In the model distribution, 96% of the 
population has a normal distribution of log P with mean value 
0.75 and standard deviation 0.18, and 4% has a normal distribu- 
tion with mean value 1 .48 and standard deviation 0.20. 



3. Modelling the Cepheids observed by Gaia 

We create a synthetic population of Galactic Cepheids with 
each object assigned a number of properties generated from the 
relevant probability distribution models. In order to study the 
Cepheid P-L relation, each object is given a period (P) and an 
absolute visual magnitude (My), as well as a V - / colour as 
required for the transformation to the Gaia wide-band G mag- 
nitude. Finally, each Cepheid is given a position in the Galaxy 
from a given distribution model. We base our modelling on the 
obse i'ved distribution of the 455 Cepheids in the lBerdnikov et all 
(120001) catalogue. 

3.1. Cepheid properties 

It is now over 100 years since the discove ry of the P-L relation 
dLeavittl [T908t iLeavitt & Pickering! [l912l) . It has since played 
an important part in the extragalactic dis tance ladder, and its 
calibration is still of great interest tod ay (ISandage et al.ll2004t 
iFouque et al1l2007 l: 'Ngeow et al]|2009l) . Because of the difficul- 
ties involved with observations in the Galactic plane and the 
larger number of Cepheids known in the Magellanic Clouds, the 
majority of the calibrations have used Cepheids in the LMC and 
SMC at least to determine the slope of the P-L relation. The cali- 
brations show a discrepancy between the P-L relations in the dif- 
ferent galaxies, which implies that the luminosities also depend 
on some other property than the pulsation period. With mea- 
sured metallicities becoming available for an increasing num- 
ber of Cepheids, recent studies have shown the metallicity to be 
a lik ely candidate, although the results are far from conclusive 
(iGro enewegen 2008; Romaniello et al. 2008). In this work, we 
generally assume all Cepheids to follow the same standard P- 
L relation, disregarding a possible metallicity effect. This gives 



the best-case scenario for the use of the Galactic Cepheids for 
the parallax verification. In Sect. I4.2l we briefly consider how a 
metallicity-dependent P-L relation would affect the results. 

Assuming a period distribution model (see below), the period 
of each Cepheid is first randomly generated and the true visual 
magnitude then follows from the assumed P-L relation. As the 
P-L relation will be calibrated simultaneously with the parallax 
zero-point, the precise P-L relation adopted for the simulations 
will affect the end r esults only minimally. We use the relation 
from lSandage et al.l (l2004 . 

= -3.087 logP- 0.914, (1) 

whe re P is given in days. The period-colour relation by 
Tammann et al ] (120031) is then used in the same way to gener- 
ate the intrinsic V - I colour, 

(V - /)() = 0.256 log P + 0.497 . (2) 

In this process we have neglected the intrinsic dispersion 
of the P-L relation due to the finite width of the instabil- 
ity strip i n the underlying period-luminosity-c olour relation 
(Madore & FreedmanI [l991l) . From LMC data (lUdalski et al.l 
1999) the dispersion is found to be about 0.16 mag in My and 
0.11 mag in M/. The Gaia G magnitude being intermediate be- 
tween V and / for typical Cepheid colours (cf. Eq. |6), the disper- 
sion in Mc is presumably interm ediate as well. When using the 
reddening-free Wesenheit index, lUdalski et al.l (1 19991) found a 
considerably smaller dispersion of 0.076 mag. For our best-case 
scenario we ignore the dispersion in the P-L relation, assuming 
(optimistically) that it can be accounted for by appropriate mod- 
elling of the full period-luminosity-colour-(metallicity)-(other 
factors) relation for the nearby Cepheids, using distances from 
Gaia. The effects of an intrinsic dispersion are similar to those 
of an uncertainty in the correction for extinction, which we do 
however investigate (Sect. 1431 1. 

For the modelling of the Cepheid periods, we use the 
iBerdnikov et al.l (l2000l) catalogue to obtain a likely period dis- 
tribution. Because the period is related to the luminosity, and in 
order to avoid observational biases, we base our model on the 71 
Cepheids within the completeness limit discussed in Sect. 13.21 
In the upper panel of Fig. |2]we plot the normalized distribution 
of periods both for the full catalogue of 455 Cepheids (shaded 
histogram) and for the volume-complete sample of 7 1 Cepheids 
(line histogram). We note that the full catalogue contains a larger 
fraction of long-period (i.e., high-luminosity) Cepheids, as can 
be expected for a sample that is at least paitially limited in ap- 
parent magnitude. The full catalogue suggests a bimodal distri- 
bution, less evident in the volume-complete sample. In order to 
reproduce both the strong short-period peak and the long-period 
tail of the distribution, we fit two overlapping Gaussian functions 
to the volume-complete sample. The resulting model period dis- 
tribution is shown as the line histogram in the lower panel of 
Fig. 12] Although the fit to the volume-complete distribution is 
far from perfect, the model reproduces the gross distribution rea- 
sonably well with a simple continuous function. 

3.2. Spatial distribution of Galactic Cepheids 

We assume the Cepheids to be distributed axisymmetrically 
around the Galactic Center, meaning that the number density is 
only a function of the galactocentric (cylindrical) radius R and 
the vertical distance z to the Galactic plane: 

N(R,z)^m)xmz), (3) 
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Fig. 3. Mean column density of Cepheids within a projected dis- 
tance r cos b from the Sun. The Berdnikov Cepheids are given in 
the solid (blue) curve, and the dashed (red) curve represents sim- 
ulated data using the radial distribution model. Note the plateau 
of constant density between 0.5 and 1 kpc indicating that the 
Berdnikov catalogue is complete within 1 kpc from the Sun. 



where £(7?) is the surface density at distance R from the axis, 
and f{R, z ) is the density d istribution perpendicular to the disk. 
Since the iBerdnikov et al.l (2000) catalogue is only complete 
within ~ 1 kpc around the Sun, it is necessary to estimate 
the radial distribution by some other means. Since classical 
Cepheids are young and massive stars, we can expect them to 
follow roughly the same distribution as other young and bright 
stars. iMcKee & Williams ( 1997) and WilHams & McKee (1997) 
found the radial distribution of OB associations to be best de- 
scribed by an exponential function Z(/?) oc exp(/?//?o) with 
^0 = 3.5 kpc, and we choose this also for the Galactic Cepheids. 
From the galactocentric radius, the position is then generated 
from X - R sin 6 and y - R cos 9, where is randomly picked 
between and 2n. 

The vertical distribution of Cepheids in the Berdnikov cat- 
alogue is found to be well fitted by a hyperbolic secant law 
dvanderKruitll 19881) . 
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f{R,z) = — sech 



(-) 



(4) 



where zo{R) is the radius-dependent scale height. From the 
Berdnikov catalogue we obta in for the solar neighbo urhood 
zo{R = 8.0 kpc) = 75 pc. lAmores & Lepind (|2005) found 
a scale height zo exp{R/12.5 kpc) for the distribution of 
Galactic HI and H2, and we therefore adopt zq{R) = (40 pc) x 
exp(7?/ 12.5 kpc) for the Galactic Cepheids. 

The surface density of Cepheids at the Sun's distance from 
the Galactic Cen tre, 2(/? = 8 kpc), can also be estimated from the 
IBerdnikov et al.l ( I2OOO) catalogue, if we assume that the Solar 
neighbourhood is representative and that all Cepheids at high 
galactic latitudes have been included in the catalogue. In Fig. [3] 
we plot the mean surface density of the Berdnikov Cepheids 
within a projected distance r cos b from the Sun (r is the distance 
from the Sun and b the Galactic latitude). We note a plateau of 
roughly constant column density ^ 2.5 x 10"^ pc"^ between 0.5 
and 1 kpc, before it falls off at larger distances. The plateau is 
believed to be real, and arises because of the relatively large ra- 
dial scale length of the Cepheid number density. The fall-off oc- 
curs at the distance where the sample is no longer complete. We 



therefore conclude that the Berdnikov catalogue is complete to a 
projected distance of about 1 kpc. 

To estimate the total number of Galactic Cepheids, we use 
the radial distribution described above and keep on generat- 
ing Cepheids until the surface density of the generated sample 
agrees with the surface density of the Berdnikov sample within 
the completeness limit. The discrepancy between the dashed 
(model) and solid (observed) curves in Fig. |3] represents all 
the Galactic Cepheids that have yet to be detected. This ex- 
trapolation results in a total of about 20000 Galactic Cepheids. 
Ch anging the radial sca l e len gth to /?() = 2.5 kpc, as presented 
by Binnev & M emfieldl (119981) for the Galactic thin disk, only 
changes the total number of Cepheids by 10% up to 22 000. 
Th ese numbers a re slightly larger than the 15 000 estimated 
by 'Maia ess et al.l |2009), who did not take the radial gradient 
into account. Our numbers are in reasonable agreement with ex- 
pectations from models of star formation and stellar evolution. 
E.g., assuming that all stars with masses above 5 Mq become 
Cepheids with a mean lifetim e of 2 Myr, a total star forma- 
tion rate of 3 Mq yr and the lSalpeteJ (Il955h IMF, we expect 
~ 20 000 Cepheids in the Galaxy. 

3.3. Gaia observations 

To simulate how the Cepheids are observed by Gaia we need 
their positions and apparent G magnitudes as seen from the 
Sun. The positions are immediately obtained from the simulated 
galactocentric (x, y, z) coordinates by subtracting the coordinates 
of the Sun. The apparent V magnitude is given by 



y = Mv + 5 log r - 5 - Av 



(5) 



wher e r is the heliocentric dist ance in pc and Ai/ the total extinc- 
tion. lAmores & Lepind (l2005h describe an axisymmetric three- 
dimensional Galactic extinction model obtained from observa- 
tions of the distribution of Galactic HI and H2 gas. We use this 
model to obtain both Ay and the colour excess in V-/, assuming 
the ra tio Ry i = Ay/Ev-i - 2.42 (Cox 2000; iTkmmann et al.l 
l2003h . The Gaia G band is similar to V for blue objects but 
brighter for red objects. The V magnitude is transformed into 
G by the following relation from Lindegren (2010): 

G = y - 0.017 - 0.088(y -/)-0.163(y- 7)2 -H0.009(y- 7)^(6) 



(cf.lJor di et al.ll20l"oh . The synthetic population is then observed 
according to a standard model developed by the Gaia Data 
Processing and Analysis Consortium (DPAC; iMignard et al.l 
2008). Gaia has a bright limit of G = 6, where the detectors sat- 
urate, and a faint limit of G = 20, corresponding to y ^ 20-25. 
The Gaia observational model predicts the standard error in the 
measured parallax, cr^, for each object in this range, taking into 
account the standard error per scan across the object, depending 
on the G magnitude, and the number and geometry of the scans 
over the five year mission, depending on the object's position on 
the sky. The observed parallax (zitg) is then obtained by adding 
a normally distributed random measurement error, with the c al- 
culated cr^, to the true parallax. 

In Fig. ID the Cepheids of the Berdnikov catalogue are com- 
pared with the synthetic Cepheid population, divided into five 
bins in apparent G magnitude. We estimate that Gaia will ob- 
serve roughly 9 000 Galactic Cepheids, or almost half of the to- 
tal population. This value is found to be relatively insensitive to 
variations in the Cepheid distribution parameters. We note that 
even though the extinction towards the Galactic Centre is be- 
lieved to be very large, several hundred Cepheids are still visible 
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Fig. 4. The Galactic Cepheids as observed from the Sun, with the Galaxy seen face on. The Sun is positioned at ix,y) = (-8,0) 
kpc, and the Galactic center at (0,0). The upper left panel shows the Berdnikov sample, and the other five panels show the synthetic 
Cepheid population divided into different G magnitude bins. The Cepheids in the lower right panel and those brighter than G = 6 
would not be observed by Gaia. 



near and even behind it. These stars are all found at |z| of hun- 
dreds of pc, where the total extinction is relatively small com- 
pared to in the plane. This can also be seen in Fig. |5] where 
the vertical distribution of Cepheids towards the Galactic Centre 
(|/| < 5°) is plotted. At distances larger than 5 kpc, there are no 
Cepheids visible in the Galactic plane. 

4. Analysis of simulated observations 

With the models described in the previous section it is possible 
to generate a list of observed Galactic Cepheids together with 
the simulated data (e.g., P, G, V - I, vjg, ctj^). In this section 
we describe the tools used to analyse the simulated data and the 
results of the parameter fitting. 

4.1. Parameter fitting 

The idea is to use the observed parallaxes of nearby Cepheids 
to determine the P-L relation, and to use this P-L relation for 
more distant Cepheids to determine the parallax bias. To avoid 
circularity, we make a simultaneous fit of the P-L relation and 
the pa rallax zero point to all the data. Following iKnapp et al.l 
(l200l the fitting is made in parallax space, where the error dis- 
tribution is symmetric around the true values. Stati stically, thi s 
is equivalent to the method of reduced parallaxes (Feasl) l2002h . 
and allows to handle correctly that some observed parallaxes are 
negative due to measurement eiTors. Such negative observed par- 
allaxes are statistically valid measurements, but cannot be con- 
verted to distances or absolute magnitudes. During the parameter 
fitting, all measured objects are therefore usable, avoiding a pos- 
sible bias due to selection (Feast & Catchpole 1997.) . Rewriting 
Eq. (|5]l and inserting the P-L relation, we get the observation 
equation: 

mc [juas] = iQS^o.iiaios P+b-v+Ay] ^ ^ ^ ^^^^^ ^ 



where a and b are the slope and zero point of the P-L relation, re- 
spectively, c represents the global parallax zero point eiTor, and 
is expected to be zero if Gaia works well. We can safely assume 
that Gaia will be able to measure P and V with negligible un- 
certainty. If we assume negligible intrinsic dispersion of the P-L 
relation and that Ay is known exactly, we have the best-case sce- 
nario for parallax zero-point verification using Cepheids. Each 
data point is then weighted entirely depending on its formal par- 
allax uncertainty, cr^. 

As a more realistic alternative, we introduce some uncer- 
tainty in the knowledge of the extinction value Ay by assuming 
a constant uncertainty cta^, - 0.05 mag for all objects. This is 
pessimistic for the bright and nearby, low-extinction Cepheids, 
but probably optimistic for high-extinction Cepheids. We then 
generate assumed extinction values that are normally distributed 
around the true values, and take the total uncertainty in Eq. Q 
to be cr'^ = [cr^ + (TO-cr/i„/2.17)^]"-^. This will lessen the im- 
portance of the nearby Cepheids, but might also introduce addi- 
tional bias effects since the calculated uncertainty uses the ob- 
served parallax and not the true one. As seen from Eq. (|7), an 
intrinsic dispersion of the absolute magnitude in the P-L relation 
will have the same effect as random errors in the assumed Ay. 
We can therefore use these experiments also to conclude on how 
such a dispersion would affect the results. 

To avoid the uncertainties associated with determining the 
extinction for each Cepheid, we also investigate the use of a 
reddening-free method equivalent to the use of the Wesenheit 
function W = V-Rv-i(V-I) ( Mado re & Freed man 1991. 20091) . 
where Ry i is the ratio of the total to selective extinction. If the 
period-colour relation is V - I - d log P + e, the observation 
equation then becomes 

where ki = a - dRy-i and k2 - b - eRy-i. It is necessary to 
introduce ki and k2 as the new unknowns since it is not possible 
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Fig. 5. The simulated inner Galaxy (|/| < 5°), with projected dis- 
tance to the Sun plotted versus height above the Galactic plane. 
The light dots correspond to Cepheids that will be observable by 
Gaia (G < 20) and the dark dots correspond to Cepheids that are 
too faint to be observed by Gaia (G > 20). 
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Fig. 6. Distribution of the fitted parameters in 1000 realisations 
of the Cepheid data. The upper three panels show the distribu- 
tions of a, b (in the P-L relation) and c (the parallax zero-point) 
for the model in Eq. (|7]l when extinction is perfectly known 
(left), and when it has an uncertainty of 0.05 mag (right). The 
bottom panels show the distribution of c for the model in Eq. (jSJ 
when Ry-i is perfectly known (left) and when it has an error of 
5% (right). The vertical Unes indicate the true values. 



to solve simultaneously for all four parameters a, b, d and e. This 
method requires however that Ry-i is known, and we investigate 
the effect of assuming a value of that is too large by 5%. 

Least-squares fitting using the Newton-Raphson iterative 
method gives the parameters a, b (or ki and k2) and c along with 
their formal uncertainties arising from the known uncertainties 
in the measured parallaxes. Biases and the total uncertainties re- 
sulting from all modelled effects can be obtained after multiple 
realisations of the Cepheid data and parameter fitting. 



4.2. Numerical results 

In Fig.|6]we present the distributions of the fitted parameters a, b 
and c after 1000 realisations of the Cepheid observational data. 
In the left panels, the correct extinction is assumed during the 
parameter fitting, meaning that the spread arises only due to the 
uncertainty in the Gaia parallaxes. This (unrealistic) best-case 
scenario leads to very well-determined P-L relation parameters 
{(Ta - 0.0035 mag dex"', cr/, = 0.0027 mag) and a parallax 
zero-point uncertainty or cr^. - 0.20 yuas. No bias is observed: 
the distributions are symmetric around the true parameter values. 
The bottom left panel shows the distribution of c when using 
the reddening-free model of Eq. ([8} with the correct value of 
Rv-i- Again, the results are unbiased but the parallax zero-point 
uncertainty is slightly larger, cr^ - 0.25 fias. 

In the right panels of Fig. |6] we show how the results are af- 
fected by an imperfect knowledge of the extinction, all other fac- 
tors being the same as in the left panels. The top three diagrams 
show the results for a, b and c when the assumed extinction has 
an uncertainty of cr^j, = 0.05 mag. The parallax zero-point un- 
certainty has increased to cr^ = 0.27 fias, and in addition there is 
a significant bias in the P-L relation zero-point (b) and a corre- 
sponding bias in the parallax zero-point of 0.3 /ias. The bottom 
right diagram shows that the reddening-free model of Eq. dH) is 
very sensitive to an error in the assumed Rv i- A 5% error in 
Rv-i introduces a bias of about 1 1 yuas in c. To keep the bias 
below 0.1 //as would require that Ry i is known to better than 
0.04%. 

The intrinsic dispersion of the P-L relation is at least about 
0. 1 mag (Sect. 13. Il l, and its expected effect on c is therefore 
at least twice as big as the 0.05 mag uncertainty in the extinc- 
tion, including a likely bias of the order of 0.6 fias. Using the 
reddening-free method reduces the scatter both due to the ex- 
tinction and the intrinsic dispersion, but instead the results are 
then very sensitive to an error in Ry-i, as we have seen. 

We have briefly inv estigated the e f fects o f a metallicity 
dependent P-L relation. iGroenewegenI (l2008b found My - 
-2.60 logP- 1.30 H-0.27[Fe/H] for the Galaxy, where the metal- 
licity dependence has an uncertainty of 0.30 mag dex"'. We as- 
sumed this metallicity dependence to be the true one, and im- 
plemented a radial metallicity gradient of [Fe/H1(j^) = 0.42 - 
0.052/? with an internal scatter of 0.1 dex (iLemasle et al.ll2008h . 
Assuming that the metallicity of each Cepheid can be deter- 
mined with an accuracy of o-[Fe/H] =0.1 dex and adding a metal- 
licity parameter in the P-L relation, the parallax zero-point un- 
certainty increases to cTc - 0.52 yuas in the case of a perfect 
knowledge of the extinction. This model also results in an in- 
creased bias of 0.45 //as in c. If metallicity is not accounted for 
in the P-L model, the scatter in c is somewhat reduced but the 
bias is even larger (1 .3 /(as). 

The experiments described above used all the Cepheids ob- 
served by Gaia in the least-squares fitting, independent of their 
individual accuracies and degrees of extinction. As we have 
seen, the results are very sensitive to extinction errors. It is pos- 
sible that this sensitivity is a consequence of including many 
Cepheids with large extinction in the analysis. In order to in- 
vestigate this we tried various ways of removing the worst data, 
e.g., by only using Cepheids within certain distance or extinc- 
tion limits. This gave only a slight improvement in terms of the 
scatter in c, but was usually found to introduce additional bi- 
ases that proved difficult to avoid. One reason for this could be 
that the selection of 'good' data depend on measured quantities, 
which in the case of noisy data invariably introduces selection 
biases. We also note that the method requires the observation of 
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both nearby and distant Cepheids in order to separate the P-L 
zero point (b) from the parallax zero point (c); excluding either 
nearby or distant Cepheids from the analysed sample introduces 
large statistical uncertainties in both parameters. 

5. Conclusions 

In order to explore the full statistical potential of the Gaia par- 
allaxes it is desirable that the global parallax zero point can be 
verified to within 0.1 yuas. We have explored the possible use of 
Galactic classical Cepheids for this purpose. 

A model of the Galactic Cepheid population has been formu- 
lated which allows us to simulate their observation by the Gaia 
satellite. From the simulated data, we have made simultaneous 
fits of the P-L relation and the Gaia parallax zero point under a 
variety of assumptions. 

We find that the parameters a and b of the P-L relation can be 
determined with a typical precision better than 0.01 mag dex"' 
orO.Ol mag, respectively, which is far better than current calibra- 
tions. The results for the parallax zero point c are less encourag- 
ing. Even under optimal circumstances (accurate knowledge of 
extinction, metallicity, etc), the Galactic Cepheid method cannot 
determine c better than to within a few tenths of a yuas. Moreover, 
we find that the resulting c is very sensitive to errors in the ex- 
tinction correction, or to an error in the Ry i value if a reddening- 
free method is used. Attempts to improve the situation, e.g., by 
limiting the sample to low-extinction Cepheids, were largely un- 
successful due to the introduction of additional biases caused by 
the selection being made from observed values. 

By extrapolating Cepheid statistics from the lBerdnikov et al] 
( 1200 0) catalogue, we estimate the total number of Galactic 
Cepheids to be ~ 20 000. We estimate that Gaia will observe 
about 9 000 of them, which is a factor ten larger than the cur- 
rently known number. Although many of them are faint, their 
observation by Gaia will greatly improve our knowledge of the 
P-L relation and its dependence of other factors such as metal- 
licity. A detailed global modelling of their characteristics is very 
worthwhile, and should take into consideration a possible par- 
allax zero point error. However, the ultimate astrophysical ver- 
ification of Gala's parallax zero-point is likely to depend on a 
combination of many different methods including the presented 
one. 

Finally, if we assume that the parallax zero point can be 
verified to a good accuracy without the use of Cepheids, one 
could use the Cepheids observed by Gaia to learn more about 
extinction. With the highly accurate observations that Gaia will 
provide, this method could yield very precise mapping of the 
Galactic extinction (R and Ac) with an accuracy that has previ- 
ously not been achievable. 
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